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The dynamics of a quantum system undergoing frequent measurements (quantum Zeno effect) 
is investigated. Using asymptotic analysis, the system is found to evolve unitarily in a proper 
subspace of the total Hilbert space. For spatial projections, the generator of the "Zeno dynamics" 
is the Hamiltonian with Dirichlet boundary conditions. 
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o : 

Frequent measurement can slow the time evolution of a quantum system, hindering transitions to states different 
■ from the initial one This phenomenon, known as the quantum Zeno effect (QZE), follows from general features 

of the Schrodinger equation that yield quadratic behavior of the survival probability at short times ||. 

However, a series of measurements does not necessarily freeze everything. On the contrary, for a projection onto 
a multi-dimensional subspace, the system can evolve away from its initial state, although it remains in the subspace 
defined by the "measurement." This continuing time evolution within the projected subspace we call quantum Zeno 
dynamics. It is often overlooked, although it is readily understandable in terms of a theorem on the QZE Q that we 
^ ■ will recall below. 

\ The aim of this article is to show that Zeno dynamics yields ordinary constraints. Under general conditions, the 
evolution of a system undergoing frequent measurements takes place in a proper subspace of the total Hilbert space 
and the wave function satisfies Dirichlet boundary conditions on the domain defined by the measurement process. 
Moreover, the evolution is generated by a self-adjoint Hamiltonian and remains reversible within the Zeno subspace. 
This shows that the irreversibility is not compulsory, as noted in Q. 

The QZE has been tested on oscillating systems || and there has been a recent observation of non-exponential 
decay (leakage through a confining potential) at short times || . Although these experiments have invigorated studies 
on this issue, they deal with one-dimensional projectors (and therefore one-dimensional Zeno subspaces): the system is 
forced to remain in its initial state. This is also true for interesting quantum optical applications tm. The present work 
(—| . therefore enters an experimentally uncharted area, although the property of being a multidimensional measurement 
is not at all exotic, and in particular applies to the most basic quantum measurement: position. The latter is the 
paradigm for the present work. 

We introduce notation. Consider a quantum system, Q, whose states belong to the Hilbert space 7i and whose 
evolution is described by the unitary operator U(t) — exp(-iHt), where H is a time-independent semi-bounded 
Hamiltonian. Let E be a projection operator that does not commute with the Hamiltonian, [E, H] ^= 0, and ETL = TLe 
the subspace defined by it. The initial density matrix po of system Q is taken to belong to TLe ■ 



X 



p =Ep E, Trp = l. (1) 
The state of Q after a series of ^-observations at times tj = jT/N (j — 1, • • • , N) is 

pW(T) = V N (T)p V^(T), V N {T)^[EU{T/N)E] N (2) 
and the probability to find the system in TLe ("survival probability") is 

P W(T)=Tr[v N (T)p vUT)}. (3) 
Our attention is focused on the limiting operator 

V(T) = lim V N (T). (4) 

N— >oo 

Misra and Sudarshan Q proved that if the limit exists, then the operators V(T) form a one-parameter semigroup, 
and the final state is 

p(T) = lim p N (T) = V(T)p V\T). (5) 

TV— >oo 
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The probability to find the system in TLe is 



V(T) = lim P {N) (T) = 1. (6) 

N—>oo 

This is the QZE. If the particle is constantly checked for whether it has remained in TIe-, it never makes a transition 
to (Wb) 1 . 

A few comments are in order. First, the final state p(T) depends on the characteristics of the model investigated and 
on the measurement performed (the specific forms of Vjv and V depend on E). Moreover, the physical mechanism that 
ensures the conservation of probabilities within the relevant subspace hinges on the short time behavior of the survival 
probability: probability leaks out of the subspace TIe like t 2 for short times. Since the infinite- N limit suppresses 
this loss, one can inquire under what circumstances V(T) actually forms a group, yielding reversible dynamics within 
the Zeno subspace. 

In this article we show that Zeno dynamics for a position measurement yields a particular kind of dynamics within 
the subspace defined by that measurement, namely unitary evolution with the restricted Hamiltonian, and with the 
domain of that (self adjoint) operator defined by Dirichlet boundary conditions. This elucidates the reversible features 
of the evolution for a wide class of physical models. As a spinoff, our proof provides a rigorous regularization of the 
example considered in 

We start with the simplest spatial projection. Q is a free particle of mass m on the real line, and the measurement 
is a determination of whether or not it is in the interval A = [0, L] C R. The Hamiltonian and the corresponding 
evolution operator are 

H = , U(t) = expHitf ) . (7) 

H is a positive-definite self-adjoint operator on L 2 (R) and U(t) is unitary. We study the evolution of the particle 
when it undergoes frequent measurements defined by the projector 



E A = I dx XA (x)\x)(x\, (8) 

where x A i s the characteristic function 

v (x) = ! 1 fo " e ^=[°^] (9) 

Xa{ > \0 otherwise ' { ) 

Thus Ea is the multiplication operator by the function \A- We study the following process. We prepare a particle 
in a state with support in A, let it evolve under the action of its Hamiltonian, perform frequent Ea measurements 
during the time interval [0, T], and study the evolution of the system within the subspace TLe a = EaH. We will show 
that the dynamics in TLe a IS governed by the evolution operator 

V(T)=exp(-iTH z )E A , (10) 

with H z =£- + V A (x), Va(x) = ( ° f " X e . A • (11) 
2m I +oo otherwise v ' 

This is the operator obtained in the limit (||). In other words, the system behaves as if it were confined in A by rigid 
walls, inducing the wave function to vanish on the boundaries x — 0, L (Dirichlet boundary conditions). 

We now prove our assertion. Let the particle be initially (t = 0) in A. We recall the propagator in the position 
representation [||,[)| 

G(x,t;y) = (x\E A U(t)E A \y) = x A {x){x\U(t)\y)x A {y) 



where t = T/N is the time when the first measurement is carried out and the particle found in Ea- To study the 
properties of G we choose a complete basis in L 2 (A) 



A ( nnx \ 

u n (x) = (x\u n ) = W — sin — J (n = l,2, ...). (13) 
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When these functions define the eigenbasis of H , H is self adjoint and 



H\u n ) = E n \u n ) , E n = -^J- , (14) 



so that H has Dirichlet boundary conditions. The matrix elements of G are 



G mn (t) = {u m \E A U(t)E A \u n ) = I dx I dyu m (x)J , exp 



2irith 



im{X ~ Vr u n (y). (15) 



2ht 



Let r = x — y, R = (x + y)/2 and A = m/2ht, so that 



G mn (X) = \l— dR dr u m {R + r/2)u n (R - r/2) exp \i\r 2 } . (16) 



where tq{R) — L — \L — 2R\. We now use the asymptotic expansion 



5(A) = \ - [ dx f(xy Xx2 = 5stat (A) + ffb ound(A), (17) 

V "T J_ a 

where 

3 s tat(A)=/(0) + ^/"(0) + O(A- 2 ), (18) 

5bound (A) = f a [/(a) + f(-a)} + 0(A- 3 / 2 ) (19) 
2iav«7rA 

arc the contributions of the stationary point x = and of the boundary, respectively. By expanding the inner integral 
in |L6|) as in (0)-(|l|) one gets 

X" r ro(R) 

dr u m (R + r/2)u n (R - r/2) exp [zAr 2 l 

= u m (R)u n (R) + [u m (R + r/2)u n (R - r/2)] r=0 + 0(A~ 3 / 2 ) . (20) 

(Note that the contribution of the boundary vanishes identically.) Using this result, we integrate by parts and after 
a straightforward calculation obtain 



G mn (t) — / dR 
Jo 



it -ft 2 d 2 

u m (R)u n (R) - T u m (R)- — -j^u n (R) 
ft 2m dR z 

„2 



+ Q(t 3 / 2 ) 



= (u m \u n ) - % -^{u m \^-\u n ) + 0{t 3 / 2 ) 
ft 2m 

l-jE^+0(t 3 ' 2 ). (21) 

With this formula we can carry out the limit required in Eq. (Q). At time T, in the energy representation, the 
propagator becomes 

(T) = (m\V(T)\n) 

= lim Y G mni (T/N)G nin2 (T/N) ■ ■ -G nN _ in (T/N) 

"l-njv-l 



This is precisely the propagator of a particle in a square well with Dirichlet boundary conditions and proves (|1 0|) ([1 ip . 
Note that the t 3 ! 2 contribution in ( pl| ) drops out in the N — > oo limit since it appears as N x 0(1/7V 3 / 2 ). It is worth 
emphasizing that although this result has been proved using the basis ( |T3| ) , the information obtained is a property of 
the propagator, and therefore holds true in general. Our choice was a matter of convenience. With a different basis 
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and nonvanishing boundary conditions, the dominant contribution of order A -1 / 2 in (/bound(A) would have given a 
nondiagonal term in (p0|)-(p2|), showing that the chosen basis is not the right eigenbasis of Hz (i.e., for the limiting 
object). 

This result can be generalized to a wide class of systems. Let 



H = -f- + V, U(t) = exp(-UH), 
2m 



(23) 



where V is a regular potential. (It may be unbounded from below, for example V(x) = Fx, although within A the 
total Hamiltonian H should be lower bounded.) The measurement performed is again application of the projector 
(H) and we study the short-time propagator 



G{x,t;y) = x A {x) 



m 



exp 



im{x — y) 2 




exp 


2ht 



it(V(x) + V(y)) 



21, 



x A (y)- 



(24) 



The basis to be used for representing the propagator is again that of the Hamiltonian with Dirichlet boundary 
conditions in [0, L] 



As before (r = x — y, R = (x + y)/2, A = m/2Ht) 



E„\u n ) , Un(x)\ x =o,L = 0. 



G mn (X) — \ — 



dR 



MR) 



-ro(fl) 



(25) 



(26) 



Using the asymptotic expansion (p"7|)-(|l9|), a calculation identical to the previous one yields 



dR 



it 

u n {R)u m (R) - —u n (R) 
n 



-h 2 d 2 

'2^dR 2 



V(R) u m (R) 



0{t^ 2 ) 



(u n \u rn ) - % ±{u n \ ( |- + V ) \u m ) + 0(t 3 / 2 ) 
a \ 2m 



= S„r 



it 

l- T E r 



+ 0{t^ 2 ) 



(27) 



and the limiting propagator at time T again reads 

G (T)=S e - lTE " /h 



(28) 

Again, the simplicity of the proof is due to the choice of the basis (^H|), satisfying Dirichlet boundary conditions. 

We have also obtained an improvement with respect to earlier approaches to this problem. The aforementioned 
theorem by Misra and Sudarshan requires that the Hamiltonian be lower bounded from the outset. However, we 
need only require that the Hamiltonian be lower bounded in the Zeno subspace. Despite the fact that for unbounded 
potentials (like V = Fx) H may not be lower bounded on the real line, the evolution in the Zeno subspace is governed 
by the Hamiltonian 



Hz= 2^ + Va{x1 Va{x) = +oo otherwise 



(29) 



that can be lower bounded in A, yielding a bona fide group for the evolution operators. 

The above calculation and conclusions can readily be generalized to higher dimensions, so long as the measurement 
projects onto a set in W 1 with a smooth boundary (except at most a finite number of points). We again take 
x, y £ M. n and let the measurement-projection be defined by A C K™, which is not necessarily bounded. Again setting 
r = x — y, R = (x + y)/2, Eq. (p6|) becomes 



G mn (A) 



i/2 



dR 



dr u n {R + r/2)e 



-itV(R+r/2)/2h 



u n (R-r/2)e 



-itV (R-r /I) jlh e i\r 2 



(30) 



D(R) 
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where D(R) is the transformed integration domain for r. The n-dimensional asymptotic expansions ([T7j)-(^9|) read 



3stat(A) = /(0) + ± A/(0) + 0(A- 2 ), (31) 



9b 



oun 



d (A) = 0(\- 1/2 ) x /(boundary) + 0(A~ 3/2 ) (32) 



and the theorem follows again because / vanishes on the boundary (Dirichlet). The proof is readily generalized 
to non-convex and/or multiply-connected projection domains, the only difficulty being that the integration domain 
in (|3^) must be broken up. It is interesting to notice that at those points at which the boundary fails to have a 
continuously turning tangent plane, the asymptotic contribution of the discontinuity in the boundary in ( |32| ) would 
be of yet higher order in A. 

In conclusion, for traditional position measurements, namely projections onto spatial regions, we have shown that 
Zeno dynamics uniquely determines the boundary conditions, and that they turn out to be of Dirichlet type. This 
is also relevant for problems related to the consistent histories approach [pd[^3[ , where different boundary conditions 
were proposed. For us, the frequent imposition of a projection, the traditional idealization of a measurement, provides 
all the decohering of interfering alternatives that is needed. On the other hand, in the works just cited one seeks 
a restricted propagator (using the path decomposition expansion fl4|| ) and such interference can occur. A second 
issue discussed in these works (especially [12Q is the validity of the Trotter product formula in certain cases. Again, 
our implicit use of this formula (in Eq. (12), etc.) is nothing more than its use for a free particle (or a particle 
in an ordinary potential). This is because the propagator of Eq. ( |l2|) provides time evolution under a sequence of 
operations: the particle evolves freely (on the entire line) for a time i, and then one applies the projection (left and 
right multiplication by the operator £U). 

The present work has implications for the notion of "hard wall," as used for example in elementary quantum 
mechanics. Everyone would agree (we expect) that this notion is an idealization. However, in many cases where this 
idealization is useful the "wall" is dynamic rather than static, the result of some fluctuating atomic presence. In 
this article we have a sufficient condition for the validity of this notion in a dynamic situation. Moreover, there is 
a quantitative framework (arising from our asymptotic analysis and finite-time-interval QZE effects) for gauging the 
effects of less than perfect hard walls. 
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